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Abstract 

We study truncated parabolic function and demonstrate that it is a representation of the Dirac 5 
function. We also show that the truncated parabolic function, used as a potential in the Schrodinger 
equation, has the same bound state spectrum, tunneling and reflection amplitudes with the Dirac 
5 potential as the width of the parabola approximates to zero. Dirac 5 potential is used to model 
dimple potentials which are utilized to increase the phase-space density of a Bose-Einstein conden- 
sate in a harmonic trap. We show that harmonic trap with a 5 function at the origin is a limit case 
of the harmonic trap with a symmetric truncated parabolic potential around the origin. Therefore, 
we propose that the truncated parabolic potential is a better candidate for modeling the dimple 
potentials. 

Keywords: Truncated Parabolic Potential, Dirac 5 function, Dimple Potential, JWKB Approxima- 
tion 
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I. INTRODUCTION 



The dimple potentials are used to increase the phase space density of the dilute Bose gases 
l|-|6j. This increase gives rise to the the formation of a Bose-Einstein condensate (BEC) 
when the initial sample is just above the critical temperature. So, the dimple potentials 
provide a means of inducing condensation without cooling [4]. Moreover, the improvements 
in the trapping techniques of cold atoms made possible to realize Bose-Einstein condensates 
in low dimensions [z-ll|. Thus, theoretical models of dimple potentials in one dimension 



have been studied in the literature 
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131 ] . While some of these models utilize an inverse 



Gaussian function others adopt Dirac 5 functions to describe dimple potentials. We will 
use a truncated parabolic potential V(x) = —Uq(1 — x 2 /a 2 ) for x < |a| to model dimple 
potentials first proposed in Ref. (3] for three dimensions. However, in this study Ma et. 
al. used a numeric simulation to explain for both the loading process and the subsequent 
evaporative cooling from a dimple potential. 

The advantage of the parabolic model compared to Gaussian model is that the 
Schrodinger equation for the truncated parabolic potential is analytically solvable, so one 
obtains the thermodynamic quantities of a weakly interacting Bose gas at least as a first 
approximation. Although the Schrodinger equation for Dirac 5 potentials are also analyti- 
cally solvable, one can not incorporate both the dimple depth and width together into these 
models. Only the multiplication of these quantities are used as the coupling coefficient of the 
Dirac 5 potential 12] . In addition to this, Dirac 5 models are only valid for tough and narrow 
dimples. However, experiments with relative wide dimples are also performed [4J. Therefore, 
the parabolic model which accommodates both dimple depth and width separately is more 
useful in describing dimple potentials. 

As a special case of point interactions, Dirac 5 potential has been taking interest for 
years. Atkinson and Crater investigated the effect of a Dirac 5 to the bound states of various 
potentials in one dimension [jjj]. Several authors made further studies on this subject 
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171 ] . Bound states of arbitrary but finite number of Dirac 6 shells in higher dimensions 
; |. The resonances of potentials with several Dirac 5 are examined by 
20| . In physical applications, Dirac 5 potential is mainly used to describe 



are also studiec 
Gadella et. al. 



19 



short range effects. For example, Erkol et. a 



interactions in atomic and nuclear physics 21 



used 5 potential for modeling short range 



221 ] . Moreover it is shown that 5 potentials 
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can also be used to find approximate solutions for tunneling and reflection amplitudes for 



23]. 



different potentials 

In this study, we suggest that the truncated parabolic function fs(x) = 3/ (4a) (1 — x 2 /a 2 ), 
for \x\ < a is a representation of the Dirac S function. Moreover, since the Schrodinger 
equation of the truncated parabolic potential is analytically solvable, we are able to show that 
the bound state spectrum, tunneling and reflection coefficients of the truncated parabolic 
potential go to those of the Dirac 5 potential as the width decreases while the depth increases. 
We also demonstrate that the spectrum of the harmonic trap with a truncated parabolic 
potential reduces to the spectrum of the harmonic potential with a Dirac 5 function in the 
same limit. So, the truncated parabolic potential model of the dimple trap includes all the 
results obtained for the non-interacting Bose gases by the Dirac 5 models for the dimple 
potential as a special case. 

The paper is organized as follows: In section [Til we present the solution of the Schrodinger 
equation for a harmonic trap with a truncated parabolic potential in one dimension and 
obtain eigenvalue equations and eigenf unctions, for the sake of completeness. In order to 
find the eigenvalues, the eigenvalue equations are solved numerically. In section ITTT1 we apply 
JWKB approximation to check the validity of numerically found eigenvalues. In section IIVI 
we use the sudden perturbation theory to find the transition amplitudes from harmonic trap 
to the harmonic trap with a dimple described by a truncated parabolic potential. In section 
fVT. we show that the truncated parabolic function provides a representation of the Dirac 5 
function. Moreover, we also demonstrate that the spectrum, the reflection and transmission 
coefficients of this potential reduces to those of the Dirac 5 potential as the width of the 
truncated parabolic potential goes to zero, when the depth times width is fixed. Finally, we 
conclude in Section IVT1 



II. HARMONIC TRAP WITH A PARABOLIC POTENTIAL 



In this section, we obtain the eigenvalues and the eigenfunctions of the Hamiltonian 



6 = £ + vii) (i) 



3 



in one dimension where the potential function is given as 



\mu 2 x 2 for |x| > a 

V{x) = < , (2) 



\moj 2 x 2 — U (l — for \x\ < 



2 

Here Uq > and a > represent the depth and the width of the truncated parabolic 
potential, respectively. 

The time independent Schrodinger equation for this potential is 

h 2 d 2 ^(x) 

-— a ^l +V (xMx)=E^(x), (3) 

where V(x) is given in the Eq. (JSJ). For Uq = 0, the Hamiltonian reduces to the harmonic 
oscillator Hamiltonian whose eigenvalues are E n = (n + l/2)fiw, n = 0, 1, 2, . . . and whose 
eigenfunctions are given in terms of the Hermite polynomials H n (y/muj /H x). However if 
Uq 7^ 0, the eigenvalues are no more (n + l/2)hu> and the eigenfunctions cannot be written in 



terms of the Hermite polynomials H n ( -i/ muo / h x), n = 0, 1, 2, . . . whose degrees are finite. 

The time independent Schrodinger equation of the Hamiltonian in Eq. ([1]) whose potential 
term is given in Eq. (|2]) is different for |x| > a and for |x| < a. For \x\ > a the Schrodinger 



equation reduces to 



h 2 d 2 ^(x) 1 2 2t , rnT , , 

+ -mw x = EW(x). (4) 



2m c?x 2 2 

Although the differential equation (jl]) is same as the time independent Schrodinger equa- 
tion for the harmonic oscillator potential, the solutions for \x\ > a differ from harmonic 
oscillator solution because of the continuity conditions on the wave functions and their 
derivatives at the points \x\ = a. 

We continue by rewriting the Eq. (j3J) in terms of dimensionless variable 



muo 
as 



1/2 



(5) 



^ + (2,-, 2 )*(,)=0 (6) 

where e = E/{huS). The general solutions of Eq.([6]) can be written in terms of parabolic 
cylinder functions [24] ]: 

dm = 2vHr ( y ( - X <1\{* z ^ + r(-i)^ (1 -_ A > /2 : 3/2; ^> 



2 r(i=i) ' _> 2'" r(=*i 



2' r(i^) v 2 y r(^ 



(7) 



where A = e — 1/2, T{x) is the well known gamma function and $(a,7;y) is the confluent 



hypergeometric function 



241 ] . The asymptotic behavior of D\(z) and D\(—z) are as follows: 



lim D\(z) = 0, lim D\(—z) = oo and lim D x (—z) = 0, lim D\(z) = oo . 

2— >00 2— >00 " 2— > — OO 2— ^— OO 

Because /^(t^) diverge at ±oo, the wave functions for x < —a (z < — (^) 1 / 2 a) and x > 
(2 > (^) 1/2 a) are 



ipi(z) = dDx(-z) 
tp3{z) = c 4 D x (z) 



(8) 
(9) 



where c± and C4 are constants that can be determined by the normalization and continuity 
conditions. 

For |x| < a, Eq. ([3]) takes the form 



h 2 d 2 m(x) 

2m dx 2 



+ 



-mu 2 x 2 — Un ( 1 



2\ 1 



V(x) = E^f(x) 



Defining 



UJ 2 + 



2U E + U 



ma 



2 ; 



£d, ^d — £d — 



1 



hwd " 2 

and changing the variable to Zd = a/hP x i Eq. ffTUj) takes a similar form with the Eq. 



(10) 



(11) 



dzj 



+ (2e d - z 2 d ) *{z d ) = 0. 



(12) 



Therefore the solutions of this equation are same as Eq. (J7|) with z replaced by z d and A 
replaced by A^. Hence, the wave functions for \x\ < a are 



ip 2 (z) = c 2 D Xd (z d ) + c 3 D Xd (-z d ) 



(13) 



where c 2 and C3 are constants which can be found by normalization and continuity conditions. 
Since the potential is an even function of x, the eigenfunctions are either even or odd (hence 
C4 = ±ci in Eq. (jHJ) and C3 = ±c 2 in Eq. (jTBl): 



CiDa(-«) 



for z < -{^Lfl 2 c 



c 2 \D Xd {z d ) ± D Xd (-z d )} for |*| < {?f)V 2 a 



(14) 



for z > (2f ) x / 2 a 
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where + and — signs stand for the even and odd eigenfunctions in the second and third 
regions of the wave functions in Eq. ([HP - Applying the continuity conditions on the wave 
functions and their derivatives, at \x\ = a, one can express c<i in terms of c\ and get equations 
which determine the eigenvalues of the Hamiltonian for the even and odd eigenfunctions, 
respectively: 

[^[D Xd (B) + D Xd (-B)]G x (A) = [G Xd (B)-G Xd (-B)]D x (A) (15) 
f^[D Xd (B)-D Xd (-B)]G x (A) = [G Xd (B) + G Xd (-B)]D x (A) (16) 

V 

for D X (A) £ and G X {A) ^ 0, 

where 

G x {z) = j-Dx{z), (17) 

and A = y^p a, B = a. 

The Eqs. f|T5|) and ffTB"]) . corresponding to the even and odd eigenfunctions respectively, 
have infinite number of roots and can be solved numerically to determine A and therefore 
the eigenvalues E x = (A + l/2)froj. 



III. JWKB APPROXIMATION 



In the previous section we obtained analytically the eigenvalue Eqs. (fl5|) . ( TT6|) and 
the eigenfunctions of the Hamiltonian in Eq. (TTJ) whose potential term is given in Eq. ()2]). 
However, to determine the eigenvalues one needs to solve the Eqs. ( !T5|) and ( IT6]) numerically 
331 ] . As Uo i.e. the depth of the potential increases, the eigenvalues of the low lying states 
decrease. As the eigenvalues decrease one needs to check the accuracy of the values obtained 
by numerical methods. For this aim we apply JWKB 34j approximation because it gives 



the exact eigenvalues 



26H29( for harmonic oscillator Hamiltonian whose potential term is 



one of the shape invariant potentials 27 1. 

The JWKB approximation method gives rise to the well known JWKB quantization 
formula to find the eigenenergies (e.g. see (29I]). 



j- 2 



x-l 



[2m {E n - V(x))} 1/2 dx=(n + ~)H 



(18) 



where X\ and x% are classical turning points found by the equality E n = V{xi j 2}. Since the 
potential function is different for \x\ < a and |x| > a the quantization formula differs for 



E n < (l/2)mu 2 a 2 and for E n > [l/2)muj 2 a 2 . We will denote the eigenenergies whose values 
are less than V(a) = (l/2)mu 2 a 2 by E n l \ For these eigenenergies, Eq. f fl8|) becomes 



J [2m {E n l) - V(x))] 1/2 dx = J J 2m E$ - [^muj 2 x 2 - U (l - ^) 



■f2 



2m 



dx = {n + -)/wr 

A, 



(19) 



where u;^ is defined in Eq. (fiTl) and = ±y {2E$ + £/ )/ ( mw J)- ^ ^ s eas y to solve Eq. 
(HHD for E { n ] : 

EM = -U + {n+ ~)hu d for n = 0, 1, . . . , n' - 1 . (20) 

Here, (n') is the number of the eigenstates whose eigenvalues are less than V(a). It is 
necessary to determine the lowest level index for the eigenstates in the upperlying region 
whose eigenvalues are larger than V(a) (See Eq. (122]) ). n' is determined by the formula 



'=L 



\mu 2 a 2 - EP 



(21) 



where the brackets denote the floor function. The JWKB quantization equation ( H~8|) for the 
energy eigenvalues of the eigenstates in the upperlying region becomes 



2m 



-muj 2 d x 2 
2 



dx + 



■t'2 



2m E, 



n(2) 



-mco 2 x 2 dx 
2 



(n + 2^ 7T ^ or n = n' ,n' + 1, . . . 



(22) 



where x 2 = \/{2E)/( muj 2 ) and n' is given by Eq. (l2Tj) . We take the integrals in Eq. (|22|) 
only for positive values of x utilizing the fact that the potential function is an even function. 
After taking the integrals the JWKB quantization equation for En > V(a) gives 



2{E n 2) + Up 



arcsm 



2(g£ 2) +!7o) 



+ mauj d \ 



2{E n 2) + Uq) 



(2) 



Ti — 2 arcsin 



2E 



(2) 



(2) 



mau^ 



2E n - 



a 2 = (n H — )fvK. 



(23) 



This equation also can be solved only numerically like Eq. ( fl5l) and ( fl6l) . However its 
solution is much easier than the Eq. ({15]) and (fl6|) in which one has to find the roots of 
infinite series. 
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TABLE I: The eigenenergies of the first few eigenstates for H = 1, m = 1/2, cj = 1, a = 3, Uq = 10. 
For these values n' = 5 and U(a) = 2.25. All the energies are in units of hio. 





The Eij 








Th p K;i crpn tit" p 


Anal vt" l c 


JWKB 


rjjfPpT-pripp 1 Analvtir- TWKR 1 




n 

u 




-8 8333 


0000 




1 

J- 


-6.5001 


-6 5000 


0.0001 




9 


-4.1681 


-4.1667 


0.0014 




Q 
O 


-1.8447 


-1.8333 


0.114 






0.4355 


5000 


0.0645 






9 5439 


9 6859 






6 


4.1711 


4.0504 


0.1207 




7 


5.2756 


5.2638 


0.0118 


F (2) 


8 


6.3845 


6.4181 


0.0336 




9 


7.5823 


7.5390 


0.0433 




10 


8.6277 


8.6381 


0.1040 




11 


9.7171 


9.7217 


0.0046 



We compare the values of the eigenenergies obtained by the JWKB approximation and 
those found by the numerical solutions of the Eqs. (jl~5]) and (|T6|) for different Uo and a 
values in the Tables [I] and HI1 In the Table [J we take h — 1, 2m = 1. In the Table [TT] (at 
the end of the manuscript), we take the true value for h and more realistic values for the 



other parameters 



30j. In these tables the first n' eigenenergies, denoted by En \ are less 
than V(a), and their JWKB values are calculated by the Eq. f )20|) and En s are larger than 
V(a) and their JWKB values are calculated by Eq. (!23l . 

The potential energy function given in Eq. ([2]) changes its form at |x| = a where its 
value is V(a). Unless the eigenenergies are close to this value, the JWKB approximation 
gives accurate results. However, as the values of the eigenenergies get close to the V(a), 
the accuracy of the JWKB approximation decreases. It is still useful to apply JWKB 
approximation in this interval because in numerical calculations one may skip some of the 
eigenstates since the difference between the successive eigenvalues are not uniform in this 
interval. The reason that the JWKB approximation gives accurate results for the ground 



S 



and first few excited eigenstates but not for the eigenstates whose eigenenergy values are 
around V(a) is the following: The eigenfunctions of the ground and first few excited states 
go to zero very rapidly as \x\ increases and they are almost zero at \x\ = a. So these 
eigenfunctions behave like they are the eigenfunctions of the shifted harmonic oscillator 
potential V(x) = —Uq + (1/2) mu^x 2 for which the JWKB approximation gives exact results. 
In other words they don't "feel" the change of the form of the potential. For the excited 
states whose eigenvalues are much larger than V(a), the JWKB approximation becomes 
again accurate. We can explain this also using the eigenfunctions of these states. These 
eigenfunctions are nonzero up to very large values compared to \x\ = a. So the region 
where the dimple potential exists \x\ < a is very narrow compared to the width of these 
eigenfunctions. Therefore, they are not much effected by the dimple potential. 



IV. TRANSITION AMPLITUDES 



In a recent paper, Garrett et. al. report the formation of a Bose- Einstein condensate 
without cooling in an anisotropic harmonic trap using dimple potentials [4]. In this study, 
they investigate the behavior of a Bose gas in a harmonic trap with "narrow" and "wide" 
dimples both for adiabatic and sudden processes. For the sudden turn of the dimple, the 
sudden perturbation theory can be used to calculate the condensate fraction. By describing 
the harmonic trap with a dimple potential by Eq. fl2]), it is possible to apply sudden pertur- 
bation theory since one can calculate the eigenfunctions for both case (only harmonic trap 
V(x) = (1/2) mu) 2 x 2 for — oo < x < oo and harmonic trap with a dimple potential). 

In this section, we will present the sudden perturbation calculations for one dimensional 
case. Before the turn on of the dim ple, the potential is only the harmonic potential and the 



eigenfunctions for this case are 28|, |29| : 

MZ) = ( -P(- 2 /2) H n {z) , (24) 

where z is given by Eq. The eigenfunctions of the potential given by Eq. ([2]) are found 
in Section [III and given in Eq. ( 1X41 . Using the sudden perturbation theory [28], we calculate 
the transition amplitudes as 



tnx(z) = \l—l Mz)^x(z)dz. (25) 



mtu 
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100 200 300 400 

FIG. 1: The change of the transition probability from the ground state of the harmonic oscillator 
Hamiltonian to the ground state of the harmonic oscillator with a dimple potential Poo with respect 
to the depth of the dimple Uq. 

In this equation one has to choose C\ in Eq. (TH1) such that *&\(z) is normalized to unity to 
find the transition probabilities as: 

PnX = \tnx\ 2 . (26) 



We present the change of the transition probabilities with respect to Uq from the ground 
state and from the second excited state of the harmonic oscillator to the ground state of 




50 100 150 200 250 

FIG. 2: The change of the transition probability from the second excited state of the harmonic 
oscillator Hamiltonian to the ground state of the harmonic oscillator with a dimple potential P02 
with respect to the depth of the dimple Uq. 
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the harmonic oscillator with a dimple for the values h = 1, m = 1/2, u = 1, a = 3 in the 
Figures [1] and El respectively. The transitions from the even eigenstates to odd eigenstates 
are zero since the wave functions of these states are even and odd functions. 

In following studies we aim to use the truncated parabolic potential to model dimple 
potentials used in Bose Einstein condensation experiments. Similar transition probability 
calculations from higher levels to the new ground state can be used for the calculation of the 
condensate fraction of a BEC in a harmonic trap with a dimple when the dimple is turned 
on suddenly 



V. DIRAC 5 POTENTIAL AS A LIMITING CASE OF THE PARABOLIC PO- 
TENTIAL 



Recent papers show that modeling dimple potentials by Dirac 5 functions reveal basic 



12 



13j . In this section we will present 



properties of a BEC in a harmonic trap with a dimple 
that the S potential is included in the truncated parabolic potential as a special case. 
We will first show that the function 



Ta ' 1 - 





for |x| < a 
for \x\ > a 



(27) 



represents the Dirac 5 function. There is an advantage of this representation compared 
to other representations in the literature (see e.g. 3JJ Appendix II). When the function 
fs(x) in Eq. (127|) is used as a potential for a Hamiltonian like in Eq. ([2]), the resulting 
time-independent Schrodinger equation is analytically solvable. 

It is clear that the representation given in Eq. (1271) satisfies the condition 5(x)dx = 1 
used in defining Dirac S function. We will show that this definition satisfies also the sampling 

h(0) for a — >■ where h(x) denotes an 



property of 6 function 



32|: rS(x)h(x)dx 



analytical function. That is 



lim 

a->0 



fs(x)h(x)dx 



lim 

a-*-0 



fs(x)h(x)dx = h(0). 



(28) 



For this aim we define the functions Hi(x), H2(x), H^{x) as 35] 
dH^x) N d 2 H 2 (x) x d 3 H 3 (x) 



dx 



h(x) 



dx 2 



h(x) 



dx 3 



h(x). 



(29) 
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We apply three times integration by parts to the expression lim a ^ /" fs(x)h(x)dx and get 



lim / f 5 {x)h(x)dx = lim i — [H 2 (a) + H 2 (-a)} - [H 3 (a) - H 3 (-a)] 



(30) 



Using the Maclaurin series expansion of H2(a), H2(—a) and H^(a), H^{—a) up to a 3 one can 
see that the right hand side of the Eq. ( 130]) is equal to h(0) in the limit a — > 0. 

It is well known that there is only one bound state of the Dirac delta potential (see e.g. 



M) 



h 2 

V s = -—crS(x). 
2m 



(31) 



The solution of the eigenvalue equation for this potential is k = a/2 where k 



-(2mE)/h. Therefore one finds [18j 



h 2 a 2 



E 



5m 



(32) 



Here the coupling coefficient of the 6 function is taken as —h 2 a/(2m) for calculational 
convenience. We will now calculate the bound state eigenvalues for the potential 



V{x) 



-U (l - 2jJ for |x| < a 



(33) 







for Ixl > a 



and present that for small a the eigenvalue of the ground state energy of this potential 
approximates to the ground state energy of Dirac 5 potential given in Eq. (|3T|) . The even 
and odd eigenfunctions of the bound states for the potential given in Eq. (1331) is: 



lp e {x) 



and 



ip Q (x) 



A e e KX 

Be [D, d { 

A e e~ KX 

A Q e KX 
Bo [D Jd ( 
-A e~ KX 



h • r ) B""i d \ 



x 



D 



id 1 



x 



for x < —a 
for |x| < a 
for x > a 

for x < —a 
for \x\ < a 
for x > a 



(34) 



(35) 



where v = ^2U /(ma 2 ), 7^ = (E + Uo)/{fw) — 1/2 and the coefficients A e , A Q , B e and B Q 
are to be determined from normalization and continuity conditions of the wave function and 
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its derivative. Applying these continuity conditions one gets the eigenvalue equation for the 
potential given in Eq. fl33|) of the even and odd eigenfunctions, respectively as: 



— = -yq^Tl) = G ^°! " G ^~^X a \ (36) 



where 7 = Ej{hv) — 1/2; -Da(-z) and Ga(^) are defined in Eqs. ([7]) and (ITT!) , respectively 
By taking into account the definition of the representation for 5 in Eq. (12?]) . the comparison 
of Eqs. ( l3~Tj) and (133]) reveals that a = 8maUo/(3h 2 ). As a — > 0, for constant t/oci, the 
eigenvalue Eq. ( 13"6"j) reduces to the Eq. (13"2"j) which gives the bound state eigenvalue of the S 
potential. One can show this by expanding the right hand side of the Eq. (1361) into a Taylor 
series. For constant Uoa, the left hand side of the equality in Eq. (136]) is order of a 3 / 4 



Therefore, we expand the right hand side in terms of ^mis/h a so that the expansion terms 
in the right hand side include terms up to a 3//4 . Doing this, we get the equality 



^TTy = - (27d+ i )V f a+ i(^). 



(38) 



Substituting v = ^/2U / (ma 2 ), ld = (E + U Q )/{hv) - 1/2 and 7 = E/(hv) - 1/2 into this 
equation we find 

^ 8mUna 2 , . 

= - 5 -jP- P9) 

Since cr = 8ma?7 /(3^ 2 ), Eq. ( |39|) is identical to the Eq. ( 132|) . which reveals that asa->0 

the representation given in Eq. f l27|) gives the same ground state eigenvalue with 5 potential 

for the solution of the Schrodinger equation. Doing a similar expansion to the Eq. ( 137]) . we 

see that this equation does not have a solution in the limit a — > 0. That means there is no 

odd bound state which proves there is only one bound state for the potential given in Eq. 

(133|) in the limit a 0. 

We also calculate the reflection and tunneling probabilities for the truncated parabolic 

potential given in Eq. ( 133]) for the scattering states. Then, we show that as a — > 0, the 

tunneling and reflection amplitudes for this potential reduce to the transition and reflection 

amplitudes for the Dirac 5 potential given in Eq. ( 13T1) . The wave function for a scattering 

state of a S potential incoming from —00 with an energy E > can be written as: 

'^ + Ree -^ {oix< _ a 

i>& s {x) ={ (40) 
T s e ikax for x > a 
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where k = \j2mE/h 2 and we take the coefficient of the wave incoming from — oo equal 
to unity for calculational convenience. The transition and reflection amplitudes for this 
potential is easily found to be 

2ik 



Rx 



2ik + a 
a 



(41) 
(42) 



2ik + a 

The wave function for a scattering state incoming from — oo with an energy E > for the 
potential in Eq. (|33|) can be written as: 



e ikax _|_ ikax f or x < _ a 

Cl D ld {jrfx) + C2 D, d (-^x) for \x\ < a 
Te lkax for x > a . 



(43) 



The norm squares of R and T given in Eq. ( 1431) give the probability of reflection and 
tunneling, respectively. These coefficients are calculated using the continuity of the wave 
function and its derivative at points x = \a\. In order to present the result for R, we define 

Fl (X, a) = -$ I, a 2 , 2 ) + 2(1 + A)$ (-± | a 2 , 2 ) (44) 

F 2 (A, a) = (k 2 + s 2 + a 2 , 4 ) $ (-± ~ a 2 , 2 ) - 2 S 2 (l + a 2 , 2 ) (1 + A)$ (~~, a 2 , 2 ) (45) 



F 3 (A,a) = 3(-l + zaA; + aV)<l> / 1 A 3 o„2„2, , , ^ A 5 . 



,-,aV + 2aV(-l + A)$ 



(46) 



2 '2' / v ' V 2 2 

F 4 (A, a) = 2a S 2 A$ (l - ~ a 2 , 2 ) + (ik + as 2 ) $ ± a 2 * 2 ) , (47) 

where $(a,7;?/) is the confluent hypergeometric function 24|. In terms of the functions 
Fi(A,a) to F 4 (A,a) the reflection amplitude is 



R 



IK 



-2iak 



[2aV(2 + A)$ (i 



2 ' 2 ' 



S 2 )F 1 (A,a) + 3$(i^,|,aV)F 2 (A, 



F 3 (A,a) F 4 (A,a) 
In order to demonstrate T, we additionally define 

F 5 (A, a) = -2aV(2 + A)$ Q - ^ |, a 2 * 2 ) $ 1 a 2 , 2 

F 6 (A, a) = $ U~ ~ aV J + 2aV(l + A)$ ( ~ aV 



(48) 



(49) 
(50) 
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In terms of the functions -F 3 (A, a) to -F 6 (A, a) 

ie~ 2iak k [F 5 (X, a) + 3$ (^, |, «V) F 6 (A, a)] 

F 3 (A,a)F 4 (A,a) " 1 j 

In the appendix, we show the change of the tunneling and reflection probabilities with respect 

to energy E, width a and depth Uq of the truncated parabolic potential. Now we will show 

that the tunneling amplitude T for the potential given in Eq. (133]) goes to T$ in Eq. (j4Tj). 

For this reason assuming Uoa is constant, we expand the numerator and the denominator of 

the right hand side of the Eq. ( IBTj) with respect to a and since we are interested in the limit 

a —7- 0, we keep only the constant terms which are independent of a. Doing this, we find 

lim T = — — . (52) 

a^o 3ik + AmU a/h 2 

Since a = 8maU / (3h 2 ), we see that T reduces to Tg given in Eq. (l4"Tj) for constant U a. 

Finally, we show that the eigenvalues of the Hamiltonian with the potential given in Eq. 

(J2J) approximates to the eigenvalues of the Hamiltonian with the potential 

1 H 2 

V(x) = -mtu 2 x 2 — - — <j§(x) (53) 
for Uq = 3h 2 a/(8ma) as a — > 0. The eigenvalues of even states of the Hamiltonian with the 



re Eq. (1331 is equal to E\ = (A + l/2)Hu) where A values are the roots of the 



potential of t 
equation 0, 3] 

V 2 ' = - (54) 

r(-|) 4' lMJ 

Here A = cr^/h/ (ma;). The odd eigenvalues are equal to the odd eigenvalues of the harmonic 
oscillator potential i.e. they are equal to E n = (n + l/2)hu for n = 1,3,5, ••• Q. As 
a 0, for constant Uoa = c, Eq. (ITBl reduces to Eq. ( 1541 and the roots A of Eq. (TIB]) go 



to A n = 2n + 1 where n — 0, 1, 2, 3, - • • . In order to show this, we first expand ^u/ud with 
respect to a and find 




3 2\ V 4 / 3 2\ -1/4 / S 2\ V 4 / 1 ,: S 2 

ma u \ I ma oj \ ( ma . or \ f 1 ma u 



l + ^_) w (^-) (i_l^_ + 0(a 3 )) (55) 



2c y v 2c / V 2c / V 4 2c 

for f/ a = c. Then we expand D\ d (B), D\ d (—B) , G\ d (B), G\ d (—B) in terms of B 



^mood/h a; G\(A) and D\(A) in terms of A = ymuj/h a in Eq. ( 1T5|) and obtain for the 
eigenvalue equation of the even eigenf unctions 

_ 8 (r(l/2)) 2 (mc) 3 / 4 3/4 /m \ ^ 1/2 r(l/2)r(-l/2) 8/4 

3 r(i^)r(i^) 2V4ft3/2 - ^2 C y w r(^)r(-|) ' 1 j 
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After rearranging the terms we get 

(57) 



2 U am 1 / 2 _ T(K X ' 

3 h?/ 2 uj 1 / 2 ~ tT x 



2> 

When we insert a = 8maU /(3h 2 ) and A = ay/H/(mu) into Eq. (1571) . we see that this 
equation reduces to the Eq. ( l54|) . Applying a similar procedure to the eigenvalue Eq. (flBjl 
of the odd eigenfunctions, we see that it reduces to 

r(-i/2)r(i/2) 



r(W(-¥: 



(58) 



as a — y for constant f/ a. In this limit A^ = —1/2, therefore the right hand side of Eq. ( 158]) 
can be equal to zero if and only if T(iy^) = oo. Since r(— n) = oo when n = 0, 1, 2, • • • , we 
get for A 

'' ~ ' = -n =^A = 2n + l for n = 0, 1, 2, • • • . (59) 



2 

Therefore, similar to the Dirac 5 case, in the limit a — > the odd eigenvalues of the potential 
given in Eq. fl2]) reduces to the odd eigenvalues of the harmonic oscillator potential with a 
Dirac 5 potential at the origin. 

The fact that the eigenvalues of the harmonic potential with a symmetric truncated 
parabolic potential around the origin reduce to the eigenvalues of the harmonic potential 
with a Dirac 5 at the origin, in the limit a — > for fixed U^a, shows that the results obtained 
modeling the d imple p otential with a Dirac 5 for a non-interacting BEC in a harmonic trap 



with a dimple 12|, ll3j is included in the parabolic model of the dimple as a special case. 



VI. CONCLUSION 



We propose, in this study, that the parabolic potential defined in Eq. (J2J) is more appro- 
priate for modeling the dimple potentials than the Dirac 5 potential used in the literature. 
Therefore, we showed that the parabolic potential includes the Dirac 5 potential as a special 
case and the wave functions, eigenvalues, tunneling and reflection coefficients of the Dirac 5 
potential can be obtained from those of the parabolic potential as a limiting case. 

In section [TT1 we have summarized the solution of the Schrodinger equation for a har- 
monic trap with a truncated parabolic potential in one dimension and obtained eigenvalue 
equations and eigenfunctions. Then, we presented the numerical solutions of the eigen- 
value equations. When the depth of the parabolic potential increases the eigenenergies of 
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the ground and first few excited states decrease to more negative values as shown in Table 
HT1 By plotting the eigenfunctions of these low lying states we have realized that as the 
eigenenergies decrease the numerical solutions become instable. Therefore, in Section IIIH 
we applied JWKB approximation to check the validity of the numerically found eigenvalues. 
The results in Tables HUTU show that the JWKB approximation and the numerical solution 
for the eigenenergies agree very well except in the transition region of the potential given in 
Eq. (J2J). So, we conclude that one can use JWKB approximation to find the eigenvalues for 
this potential when the numerical solutions fail. 

In section IIV1 we obtained the formula for the transition amplitudes and found the tran- 
sition probabilities from the ground and second excited eigenstates of the harmonic trap to 
the ground state of the harmonic trap with a dimple described by the truncated parabolic 
potential for different dimple depths, using the sudden perturbation theory. In the following 
studies we aim to use these and similar results in three dimensions for modeling the sudden 



turn on of dimple potential experiments {4]. 

In section |V] we have first shown that the truncated parabolic function provides a rep- 
resentation of the Dirac 5 function. Moreover, we have also demonstrated that the bound 
state spectrums of the potentials given in Eqs. fl33l) and ([2]) reduce to the spectrum of the 
Dirac 5 potential and harmonic potential with a Dirac 5 as a — > 0, for fixed Uqcl. As one can 
see from the Bose-Einstein distribution (rij) = \/{e^ Ei ~^ — 1) for a non-interacting Bose 
gas the effect of the trapping potential to the thermodynamic properties of the gas comes 
only through the eigenvalues of the potential. So, we come to the following conclusion: The 
fact that the eigenvalues of the truncated parabolic potential reduce to the eigenvalues of 
the Dirac 5 potential in the limit a — > 0, for fixed U^a, shows that the results obtained mod- 
eling the dimple potential with a Dirac 5 for a non-interacting Bose-Einstein condensate in 



a harmonic trap with a dimple 12|, [l3| is included in the parabolic model of the dimple as 
a special case. 
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|R| 2 , |T| 2 




E 



FIG. 3: The change of the tunneling and reflection probabilities |T| 2 and \R\ 2 for the truncated 
parabolic potential, Eq. (|33p . with increasing energy of the incoming wave and for fixed a = 3, 
Uq = 10 in natural units. The solid and dashed lines represent \T\ 2 and |i?| 2 , respectively. 

VII. APPENDIX 

We present in this appendix, the change of the tunneling and reflection probabilities with 
respect to energy E of the incoming wave, width a and depth Uq of the potential given in 
Eq. (I33p . The scattering wave function, reflection (R) and tunneling (T) amplitudes for 
this potential is given in Eqs. flSJ), (JI8]) and (loTj) . respectively. We first present the change 
of the tunneling and reflection probabilities with increasing energy E of the incoming wave 
in Fig. [3j for fixed potential width a = 3 and depth Uq = 10. Here we use natural units 
h — 1, to = 1/2. One can see from the Fig. [3J as E increases from zero, the tunneling 
probability first increases almost up to certainty. However, then it begins to oscillate with 
increasing amplitudes and widths such that local maximums are unity. Although the local 
minimums of these oscillations decrease almost up to zero, they never become equal to zero 
up to E = 1000. 

Then, in Fig. HI we show the variation of the tunneling and reflection amplitudes for 
increasing a and fixed Uq = 10, E = 1. We see that the tunneling probability oscillates 
such that local maximums touch one and as a increases the amplitude of the oscillations 
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|R| 2 , |T| 2 




FIG. 4: The change of |T| 2 and |i?| 2 for the truncated parabolic potential, Eq. (|33p . with increasing 
potential width a and for fixed E = 1, Uq = 10 in natural units. The solid and dashed lines represent 
\T\ 2 and \R\ 2 , respectively. 



|R| 2 , |T| 2 




FIG. 5: The change of |T| 2 and |i?| 2 for the truncated parabolic potential, Eq. (133 j) . with increasing 
potential depth Uq and for fixed E = 1, a = 3 in natural units. The solid and dashed lines represent 
|T| 2 and |i?| 2 , respectively. 

decrease. Finally, we demonstrate the change of the |T| 2 and |i?| 2 for increasing but for 
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fixed E = 1 and a = 3 in Fig. |5j Similar to the case of varying E, with increasing Uq, \T 
oscillates with increasing amplitudes where local maximums are again one. 
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TABLE II: The eigenenergies of low lying eigenstates for to = 27r20Hz, m = 23 amu, a = 11/um, 
Uo = 1.0 10~ 30 J and h is in SI units. For these values n' = 15 and V(a) = 2.75 ftw. All the energies 
are in units of huj. 



The Eigenenergies 




The Eigenstate No 


Analytic 


JWKB 


Difference | (Analytic- JWKB) | 







-72.7948 


-72.7948 


0.0000 




1 


-67.4650 


-67.4650 


0.0000 




2 


-62.1353 


-62.1353 


0.0000 




3 


-56.8055 


-56.8055 


0.0000 




4 


-51.4758 


-51.4758 


0.0000 




5 


-46.1461 


-46.1461 


0.0000 




6 


-40.8163 


-40.8163 


0.0000 


ci(i) 


7 


-35.4866 


-35.4866 


0.0000 




8 


-30.1570 


-30.1569 


0.0001 




9 


-24.8278 


-24.8271 


0.0007 




10 


-19.5005 


-19.4974 


0.0031 




11 


-14.1807 


-14.1676 


0.0131 




12 


-8.8874 


-8.8379 


0.0495 




13 


-3.6816 


-3.5082 


0.1734 




14 


1.2167 


1.8216 


0.6049 




15 


4.8125 


4.4790 


0.3335 




16 


6.0940 


6.0210 


0.0730 




17 


7.3188 


7.4137 


0.0949 




18 


8.8970 


8.7300 


0.1670 




19 


9.8991 


9.9980 


0.0989 




20 


11.3352 


11.2316 


0.1036 




21 


12.4303 


12.4396 


0.0093 




22 


13.6232 


13.6273 


0.0041 














499 


497.1836 


497.1835 


0.0001 




500 


498.1856 


498.1857 


0.0001 
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